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Abstract 
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1 Introduction 



Stochastic calculus is the most widely used mathematical formalism for mod- 
eling financial instruments , followed by the use of partial differential equa- 
tions [[| . The Feynman path integral is a formalism based on functional inte- 
gration and is widely used in theoretical physics to model quantum (random) 
phenomenon ||; it is also ideally suited for studying stochastic processes aris- 
ing in finance. In || techniques from physics were applied to the study of 
finance and in [Q the problem of the pricing of stock options with stochastic 
volatility was studied using the formalism of path integration. 

In this paper, the path integral approach is continued into the field of 
interest rates embodied in the modeling of Treasury bonds. The complexity 
of this problem is far greater than that encountered in the study of stocks 
and their derivatives; the reason being that a stock at a given instant in time 
is described by only one stochastic variable undergoing random evolution 
whereas in the case of the interest rates it is the entire yield curve which is 
randomly evolving and requires infinitely many independent variables for its 
description. The theory of quantum fields || has been developed precisely to 
study problems involving infinitely many variables and so we are naturally 
led to the techniques of quantum field theory in the study of the interest 
yield curve. 

Treating all the forward rates as independent random variables has also 
been studied in [0, || |9|, [L0| using the formalism of stochastic calculus. In 
this approach a stochastic partial differential equation in infinitely many 
variables is written. The approach based on quantum field theory is in some 
sense complimentary to the approach based on stochastic partial differential 
equations since the expressions for all financial instruments are formally given 
as a functional integral. One advantage of the approach based on quantum 
field theory is that the introduction of non-linearities as well as stochastic 
volatility is easily incorporated. 

The HJM-model |TT|is taken as the starting point of this paper. In Section 
2 the HJM-model is re-expressed in terms of a path integral, and the condition 
of no arbitrage is re-derived in this formalism. To make the formalism more 
transparent and accessible to readers not familiar with path integration, the 
well-known results for the price of futures of zero-coupon bonds as well as the 
price of a European call option and a cap for a zero-coupon bond is derived 
in Sections 3 and 4 respectively. Another more important reason for these 



1 



re-derivations is that the prices of these derivatives are expressed in a form 
which can be directly generalized to the case when we model the evolution 
of the forward rates using quantum field theory. 

In Section 5, the HJM-model is generalized to the case with independent 
fluctuations of all the forward rates; the theory is then seen to consist of 
a free (Gaussian) two-dimensional quantum field theory. The generalized 
model has a new parameter which determines how strongly it deviates from 
the HJM-model. The condition of no arbitrage is derived for the generalized 
model. 

In Section 6, the formulae for the prices of futures and options of zero- 
coupon bonds are obtained explicitly for the Gaussian quantum field theory. 

In Section 7 some conclusions are discussed as well as possible future 
directions of research. 



2 Path Integral Formulation of the HJM-model 

Bonds are financial instruments of debt which are issued by governments and 



corporations to raise money from the capital markets [12, 13]. Bonds have 
a predetermined (deterministic) cash flow; a Treasury bond is an instrument 
for which there is no risk of default in receiving the payments, whereas for 
corporate bonds there is in principle such a risk. A Treasury zero- coupon bond 
is a risk-free financial instrument which has a single cash-flow consisting of 
a fixed pay-off of say $1 at some future time T; its price at time t < T is 
denoted by P(t,T), with P(T,T) = 1. 

A Treasury coupon bond B(t, T) has a series of predetermined cash-flows 
which consists of coupons worth c, paid out at increasing times Tj's, and with 
the principal worth L being paid at time T. B(t, T) is given in terms of the 



zero-coupon bonds by |12 



K 



B(t, T) = £ <kP{t, Ti) + LP(t, T) (2.1) 

From above we see that a coupon bond is equivalent to a portfolio of 
zero-coupon bonds. Hence, if we model the behaviour of zero-coupon bonds, 
we automatically have a model for coupon bonds as well 

Consider the forward rate f(t,x), which stands for the spot (overnight) 
interest rate at future time x for a contract entered into at time t < x. The 
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price of a zero-coupon bond with the value of $1 at maturity is given by 



P{t,T) = exp{- j dxf{t,x)} (2.2) 

Note from its definition, the spot rate for an overnight loan at some 
(future) time t is r(t) and is given by 

r(t) = f(t,t) (2.3) 

The forward rate is a stochastic variable. In the K-factor HJM-model 
O , [12| , p~3f the time evolution for the forward rates is given by (sum over all 
repeated index) 

df 

-£(t,x) = a(t,x) + a i (t,x)W i {t) (2.4) 
at 

where a(t,x) is the drift velocity term and <Ti(t,x) is the deterministic 
volatility for the forward rates. From eqn.( |2.4j ) we have 



f(t, x) = /(t , x) + f dt'a(t', x) + f dt'aiit', x)W i (t') (2.5) 

J to J to 

The initial forward rate f(to,x) is determined from the market, and so 
are the volatility functions <Ji(t,x). 

Each stochastic variable Wi(t),i = 1,2. ..if is an independent Gaussian 
white noise given by 

E{W i (t)W {t')) = 5 l ,5{t-t l ) (2.6) 

Note that the forward rates f(t, x) are driven by random variables Wi(t) 
which gives the same random 'shock to all the forward rates; the volatility 
function a(t, x) weighs this 'shock' differently for each time t and each x. 
It is precisely this feature which we will generalize later such that f(t,x) is 
taken to be an independent random variable for each x and each t. 

To write the probability measure for Wi(t) we discretize t = me, with 
m — 1, 2....M = [-], and where t takes values in a finite interval depending 
on the problem of interest; then the probability measure is given by 



M K 

r[w] = n n e " fwf(m) ( 2 - 7 ) 

m=l «=1 
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(t*,tO) 



t:Time W(t) 

(tO,tO) 

x:Time to Maturity 

Figure 1: Domain of W(t) 

M K 



dW = nri\fc/ dWi(m) (2.8) 

J m=li=l V Z7V J -°° 



For notational simplicity we take the limit of e — > 0; note that for pur- 
poses of rigor, the continuum notation is simply a short-hand for taking the 
continuum limit of the discrete multiple integrals given above. We have, for 
ti <t <t 2 



V[WMM - e So (2.9) 

S Q = S Q [W,h,t 2 ] = f 2 dtW^W^t) (2.10) 

2 i=1 J tl 

dW -> j DW (2.11) 

The 'action' functional So is ultra-local with all the variables being decou- 
pled; generically, / DW stands for the (path) integration over all the random 
variables W(t) which appear in the problem. The integration variables W(t) 
are shown in Fig.(^). 

A path integral approach to the HJM- model has been discussed in 
the action they derive is different than the one given above since they use 
a different set of variables and end up with an action involving the time 
derivatives of their variables. 

A formula that we will repeatedly need is the generating functional for 
W given by the path integral 
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(t*,T) 



(tO,T) 



x:Time to Maturity 
Figure 2: Trapezoidal Domain T 



= e hS£<*Mt)m / 212 ) 

We now derive the no arbitrage condition on the drift velocity a(t, x). 
Recall having a martingale measure is equivalent to having no arbitrage for 
the price of the coupon and zero-coupon bonds ||15|| . The martingale condition 
states the following: suppose a zero-coupon bond which matures at time T 
has a price of P(t*, T) at time t* and at time to < t* has a price P(to, T); then 
the price of the bond at t*, evolved backward to time t and continuously 
discounted by the risk- free spot rate r(t) must be equal to the price of the 
bond at time to- 
rn other words, the martingale condition on the zero-coupon bond using 



(2.5) is given by 



P(t ,T) = E t0 [e J *o rm P(U,T)} (2.13) 

= P(t ,T)e-/r J DWe~fr a ^ w ^e So ^ (2.14) 

where the trapezoidal domain T is given in Fig.(|]) and 

/ = f* alt [ T olx (2.15) 
JT Jto Jt 

We can set t* = T in above equation; this will change the domain to 

a (right isosceles) triangular domain A given in Fig.(|3|) and is the largest 
domain in the problem. We have, using P(T,T) = 1, the following 



5 



t: Time 



(tO,tO) 




(T,T) 



(tO,T) 



x:Time to Maturity 
Figure 3: Maximum Domain A 



P(t ,T) = E t0 [e'^ rm ] 

e / A «(*.*) = f DWe -J A ^)w(t) e s [w] 



On performing the W-integrations we obtain from (|2.14| ) 



(2.16) 
(2.17) 



e J r 



i [** dt[f T dxa(t,x)] 2 



Dropping the integration over t we obtain [12 



(2.18) 



dxa(t, x) 



K 

Y\I dxai(t,x)f 
2 l=1 Jt 



or equivalently 



a(t,x) = (7i(t,x) J dyai(t,y) 

: Condition for No Arbitrage 



(2.19) 



(2.20) 



We have see that, as expected, the martingale condition leads to the 
well-known no arbitrage condition on the drift velocity of the forward rates. 
Consider the 2-Factor HJM-model with volatilities given by 



ai(t,x) = a 1 ;a 2 (t,x) = a 2 e x[x l) 



(2.21) 
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The no arbitrage condition given in eqn.( |2.20| ) for this case yields 



a(t, x) = a\{x -t) + ^e- x ^ (1 - e~ x ( x ^) (2.22) 

A 

3 Futures Pricing in the HJM-Model 

The future and forward contracts on a zero-coupon coupon bond are instru- 
ments that are traded in the capital markets [|12|, [yj. The forward and future 
price of P(t, T), namely F(to, £*, T) and Tito, t*, T) respectively, is the price 
fixed at time to < t* for having a zero-coupon bond delivered to the buyer at 
time t*. 

The difference in the two instruments is that for a forward contract there 
is only a single cash flow at the expiry date of the contract t*. For a futures 
contract on the other hand there is a continuous cash flow from time to to 
t* such that all variations in the price of P(t + dt, T) away from P(t, T), for 
to < t < is settled continuously between the buyer and the seller, with a 
final payment of P(t*,T) at time t*Jl^, 0. If the time-evolution of P(t,T) 



was deterministic, it is easy to see that the forward and futures price would 
be equal. 

It can be shown that the price of the futures T is given by |12| 



F(t ,U,T) = E t0 [P{U,T)\ (3.1) 
From eqns.( P3| ) and ( |2.9[ ) we have 

F{t ,U,T) = [ DWe-£ dxf(t *' x) V[W,t ,U} (3.2) 

= F(t Q ,U,T) emptor (3.3) 
where the forward price for the same contract is given by 

F(M.,T) = Sr^ (3-4) 

The trapezoidal domain T splits into a triangle and a rectangle shown in 
Fig. (H) and yields 

T = A ®n (3.5) 
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Figure 4: Domain T=Domain Aq© Domain 1Z 



The futures price is defined over the rectangular domain 71 given in Fig. (|) 
and 



We have 



with 



'* dt [ T dx (3.6) 
n J t Jt, 



expft^ = e n e-k a{t ' x) (3.7) 

J DWe-h°^ w ^e s ° (3.8) 

ex p i\ E f d A f T dx <^ x )f) ( 3 - 9 ) 

2, i=1 Jt Jt, 



where eqn. fl3.9|) has been obtained by performing the path integration 
over the W- variables using eqn. ( |2.12| ). 

Collecting results from above we have, using the no arbitrage condition 
from eqn. ( J2.2U ) and after some simplifications 



/■t* rt, rT 

flj?(to,t*,T) — — I dt dx<Ti(t,x) / dx'<Ti(t,x') (3.10) 

J to Jt Jt, 

As is expected, the future and forward prices of the zero-coupon bond are 
equal if the volatility is zero, that is, the evolution of the zero-coupon bond 
is deterministic. 

Consider the 2- Factor HJM- model with volatilities given in eqn. (|2.21|) . 
Equation ( |3.1U| ) yields 
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ttr(to,U,T) = — a x (T — t*)(t* — t ) 

- ^(l-e-^-^Xl-e-^*-^) 2 (3.11) 

which is the result given in JT|. 

4 Option and Cap Pricing in the HJM-Model 

Suppose we need the price at time to of a derivative instrument of a zero- 



coupon bond P(t, T) for a contract which expires at t* [ 12] , [13] . For concrete- 
ness we study the price of a European call option on a zero-coupon bond, 
namely C(to, t*, T, K); the option has a strike price of K and exercise time 
at t* ^> t(). 

The final value of the option at to = t* i s > as required by the contract, 
given by 



C{U,U,T,K) = (P(t*,T)-K)+ (4.1) 

= (P(U,T) - K)9(P(t*,T) - K) (4.2) 

where the step function is defined by 

{1 for x>0 

| forx=0 (4.3) 
for x<0 
For t < we have the price of C given by 

C(t , U ,T,K) = E t0 [e~ -C dtm) (P(U, T) - K) + \ (4.4) 



The expectation value in eqn. (|4.4! ) is taken by evolving the payoff function 
(P(t*,T) — K) + backward to to, discounted by stochastic spot rate r(t) = 

Using the identity 

1 r+oo 

5{z) = — / dpe*» (4.5) 
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we can rewrite eqn. (|4.4 ) as 



/+oo 
dG^(G,t*,T)(e G - K)+ (4.6) 
-oo 



where 



- //* dtf(t,t) 



*(G,U,T) = E t0 [e J *o 5(P(U,T)-e u )\ (4.7) 

= f^eV^l (4.8) 

A = lnF(t ,U,T) (4.9) 

Using (|2.2|)and (|2.5| ), we have the following 



f DWp'^o ^W^+'p k^) w ^) p S 

Note the interplay of the subdomains A and 1Z in determining the price 
of the option. Using ( ^.12| ) to perform the integrations over W yields, after 
considerable simplifications and using the no arbitrage condition given by 
eqn. ( |2.19| ) the following 

A = -y(p 2 + ^) (4.H) 



<Z 2 = E I** dt[[ T dxai(t,x)} 2 (4.12) 

i=1 Jto JU 



with 

/ dt[ 

Ho H 

To obtain eqn.( |4.1lD we have used the identity 

/ dt[ ( dxa(t,x) — [ dx<Ji(t,x) [ dy<Ji(t,y)] = -q 2 (4-13) 
Jto Jt* Jt Jt, 2 

Performing the Gaussian integration in eqn.( [4.8|) yields 
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X 



(4.14) 



Hence from above and (|4.6|) we recover the well-known result fli~7 , 18 1 that 
the European option on a zero-coupon has a Black-Scholes like formula with 
volatility given by q. 

For the 2-Factor HJM- model given in eqn.(|2.21|) we have 



q 2 = at(T-UY(U-t ) 

+ ^(l- e -Mr-t.))2(i_ e -2A(t.-to)) ( 4 . 15 ) 

A cap is a financial instrument for reducing ones exposure to interest rate 
fluctuations, and guarantees a maximum interest rate for borrowings over a 
fixed time. A cap fixes the maximum interest for a fixed period from £* to 
t* + T; the buyer of the instrument then pays for this period the maximum of 
the Libor (London interbank offered rate) L = L(i*, £* + T) or the cap rate 
K. The cap is exercised at time and the payments are made, in arrears, 
at time £* + T. Let the principal amount be V; the value of the cap at time 
i* is then given by 

Cap{U, U, T) = VT{L - K)+ (4.16) 
We have in terms of the forward rates [HI 



TL(U,U + T) = eC* -l (4.17) 
= — — - 1 (4.18) 

The price of the cap at time t < t* is given by 

Cap(t ,U,T,X) = VE t0 [e -C dtm) ( — i— - 1 - TK) + ) 

(4.19) 

/+oo 
dHT(H) (e H - 1 - TK) + (4.20) 
-oo 
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Carrying out an analysis similar to the one done for the pricing of the 
European call option we obtain, as in ( f4.14|) (note the minus sign of H) 



T(H) 



(4.21) 

J t T+t *dxf(t ,x)-?f ) 2 (4.22) 



1 1 



(-H + 



\ 2 ^cap 2 Qc< 



cap 



with q cap for the Two-Factor model given similar to (^4.15|) by 



cap 




(4.23) 



The formula above shows that a cap is equivalent to a European put 
option on the zero-coupon bond P(t*,t* + T). For a caplet the time T is 
taken to be small so that L(t*, + T) ~ /(£*, and in eqn. (|42D we have 



5 Quantum Field Theory of Treasury Bonds 

As mentioned earlier, in the HJM-model the fluctuations in the forward rates 
at a given time t are given by 'shocks' which are delivered to the whole curve 
f(t,x) by a single random variable W(t); clearly, a more general evolution 
of the instantaneous forward rate would be to let the whole curve evolve 
randomly, that is let all the forward rates fluctuate independently. The only 
constraint on the random evolution is that for every instant there be no 
arbitrage for the forward rates. 

Since there are infinitely many forward rates, we need an infinite num- 
ber of independent variables to describe the random evolution of the yield 
curve. As mentioned earlier the generic quantity describing such a system 




12 



is a quantum field 0; for modeling Treasury Bonds we will need to study a 
two-dimensional quantum field in a finite Euclidean domain. 

We consider the forward rates to be a quantum field] f(t,x) is taken to 
be an independent random variable for each x and each t. For notational 
simplicity we keep both t and x continuous; in Appendix A, the case with 
both t and x discrete is analyzed and the continuum limit discussed in some 
detail. 

For the sake of concreteness, consider the price of an options contract 
(at time to) of a T-maturity zero-coupon bond ; let the contract expire at 
time £*. Since the all the forward rates are for the future, we always have 
x > t; hence the quantum field f(t,x) is defined on the trapezoidal domain 
T bounded by t = to, t = x, t = i* and x = T, as shown in Fig (g). 

We introduce a new parameter \i which quantifies the strength of the 
fluctuations in the time-to-maturity direction x; we expect that in the limit 
of /x — > 0, we should recover the HJM- model. The simplest term which 
can control the fluctuations in the x-direction is the gradient of f{t,x) with 
respect to x. The action given in eqn.( |2.lOD is generalized to 



S[f] = C[f) (5.1) 



r 



r lf] 1 rr ^-^h 2 , ifgf ¥-«M u' 

[Jl " 2(T-t)Ll a (t,x) J fi^dx^ a(t,x) n 



(5.2) 



The initial condition is given by 



t = to, to < x < T : f(to, x) : specified (5.3) 

and the field values of f(t,x) on the rest of the boundary points of the 
trapezoid T are arbitrary and are integration variables. The quantum field 
theory is defined by integrating over all configurations of f(t, x) and yields 

Z = ( Dfe s[f] (5.4) 



Df = n / w.*) ( 5 - 5 ) 

fj. \ T J —OO 

(t,x)eT 
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Note that e s ^/ Z is the probability for different field configurations to 
occur when the functional integral over f(t,x) is performed. 

The presence of the second term in the action given in ( |5.1| ) seems to 
be justified from the phenomenology of the forward rates [|19| and is not 
ruled out by no arbitrage . This term in the action also implies that all the 
random configurations of f(t,x) which appear in the path integral ( |5.4|) are 
continuous functions of x. Forward rates that are usually observed in the 
market are continuous 0. 

However the random configurations for the forward rates are nowhere 
differentiable. It is noted in J7J that there is no evidence to indicate whether 
the actual forward rates are differentiable or not. 

The action given above is suitable for studying formal properties of the 
forward rates. However it is often simpler for computational purposes to 
change variables. Let A(t,x) be a two-dimensional quantum field; we use 
the HJM-change of variables to express A(t, x) in terms of the forward rates 
f(t, x), namely 

df 

-^-(t,x) = a(t,x) + a(t,x)A(t,x) (5.6) 
at 

The Jacobian of the above transformation is a constant and hence can be 
ignored. The action can now be written as 



= J T m (5.8) 

The quantum field theory is defined by a functional integral over all vari- 
ables A(t, x); the values of A(t, x) on the boundary of T are arbitrary and 
are integration variables; this yields the partition function 

Z = J DAe s[A] (5.9) 

Note that eqns. ( |5.6D and ( |5T7D can easily be generalized to the K-Factor 
case. It is shown in eqn.( [A.l"5|) that if we define 

W{t) = ^—J^ dxA{t,x) (5.10) 
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then for /x — > we have 

S[A) 
DA 



s =~i r dtw\t) 
i Jt 



DW 



(5.11) 
(5.12) 



From eqn.( [2.10D and above we see that we recover the HJM- model in 
the — > limit. We see from eqn.( |5.10| ) that the HJM-model is a drastic 
truncation of the full field theory and only considers the fluctuations of the 
average value of the quantum field A(t, x); it in effect 'freezes-out' all the 
other fluctuations of A(t, x). 

If one thinks of the field A(to, x) at some instant to as giving the position 
of a 'string' |7], [8], then in the HJM-model this string is taken to be a rigid 



string. The action S[A] given in (|5.7| ) allows all the degrees of freedom of the 
field A(t ,x) to fluctuate independently and can be thought of as a 'string' 
with string tension equal to -\\ in this language the HJM-model considers 
the interest yield curve to be a string with infinite tension and hence rigid. 

The moment generating functional for the quantum field theory is given 
by the Feynman path integral as 



Z[J] = - [ DAeR * e s[A] (513) 
Z J 

We evaluate Z[J] exactly in Appendix B, and from eqn.( p,17|) 



by 



Z[J] = exp - f U dt [ T dxdx'J(t, x)D(x, x'\ t, T)J(t, x') (5.14) 
2 Jt Jt 

where the propagator D(x, x'\ t, T) is given from eqn. (p.l9| ), for (3 = T — t, 



D(x, x'\ t, T) 



sinh d (/i/3) 



sinh jj,(T — x) sinh jj,{x' — t){l + sinh 2 (/i/3)6 , (x — x')} 



+ sinh /i(T — x') sinh jj,{x — t){l + sinh 2 (fj,(3)6(x f — x)} 
+ cosh(/i/?){sinh /i(x — t) sinh fi(x' — t) 
+ sinh n(T — x) sinh /x(T — x')} 



(5.15) 
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To understand the significance of the propagator D(x,x';t,T) note that 
the correlator of the field A(t, x), for to < t,t' < t* < T, is given by 



E(A(t,x)A(t',x')) = - I DAe s[A] A{t,x)A{t',x') (5.16) 
= 5(t-t')D(x,x';t,T) (5.17) 

In other words, D(x, x'; t, T) is a measure of the effect of a value of field 
A(t, x) at maturity x on its value at another maturity x' . 

Since D(x,x';t,T) looks fairly complicated, we examine it in a few ex- 
treme limits. In the limit of \i — > we have 

D(x,x';t,T) = l + 0(/x 2 ) (5.18) 

We see that, as expected, all the fluctuations in the x direction are exactly 
correlated; in other words the values of A(t, x) for different maturities are all 
the same. Defining 

]{t) = / dxJ{t,x) (5.19) 



we have from eqns. (|5.14j) and ( |5.18|) that 



Z[j] = expi f* dtf{t) (5.20) 
2 Jto 

which is the result obtained earlier in eqn. (|2.12 ). 
For /i^oowe have 

D(x, x'; t, T) ~ i^e-" 1 *-*' 1 (5.21) 

The propagator above has a simple interpretation; if the field A(t, x) has 
some value at point x, then the field at 'distances' x — < x' < x + yT l 
will tend to have the same value, whereas for other values of x' the field 
will have arbitrary values. Hence we see in this limit that the fluctuations 
in the time-to-maturity x direction are strongly correlated within maturity 
time which is the correlation time of the forward rates. 

We now derive the no-arbitrage condition for the action S[A\. Eqn. fl2.r3D 
for the martingale is unchanged; generalizing eqns. ( [2.14 ) and fl2.18| ) we have 
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exp|a(t,x) = ^ J DAe-Sr^ A ^ e Sr^ (5.22) 

- r dt r 

2 Jt Jt 



1 /■** /" T 

exp - / cit / dxdx'a(t,x)D(x,x';t,T)a(t,x') 

2 A, 



(5.23) 



Hence we have 



dxa(t, x) = - l dxdx'a(t, x)D(x, x';t,T)a(t, x') (5.24) 



n 2 Jt 

which is the generalization of eqn. fl2.19|) , and that of eqn. ( |2.20|) is given 

by 



a(t,T) = a(t,T) J dx'D(T,x';t,T)a(t,x') 

1 f T , , / / ,dD(x,x';t,T) , 
+ -J dxdx'a{t,x) V dT J a(t,x') (5.25) 



From the empirical study of forward rate curves, there is evidence [|19j that 
the naive HJM-model no arbitrage for the drift term a(t, x) is not adequate 
since it is quadratic in the volatility; in |TI| an additional term is added 



which reflects the market price of risk. In the approach of field theory, the 
additional term involving the derivative of the propagator could provide a 
better model of no arbitrage for the drift term. 
For /i^oowe have 

a(t,T) = -(T -t)a 2 (t,T) + - f dxa 2 (t,x) (5.26) 
2 2 Jt 

Note the expression for a(t, x) given above is quite dissimilar from that of 

the HJM-model given in eqn. (|2.20| ), which is the case for /i = 0; the values of 

a(t, x) given in Q5.25p for /i ^ continuously interpolate between the extreme 

values of jjt = and ji = oo. 



For the Two- Factor case given in eqn. (|2.21|) , we can exactly solve for 
a(t, x) in terms of the volatilities; the expressions are long and cumbersome. 
For the case of the One-Factor model with only a± ^ 0, we have the exact 
result that 
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a{t,x) = of(x-t) (5.27) 

which is independent of fi and the same as the HJM-model. This result 
can be seen directly from the functional integral; since 0\ =constant, in the no 
arbitrage eqn. ( |5.22j ) for ait, x\ we see that it only couples to dxAit, x) = 
(T — t)W(i); a change of variables then shows that a(t,x) does not couple to 
(a, and hence the simple result. As we will see later, the One- Factor model has 
non-trivial dependence on \x for other quantities such as futures and options. 

We have from eqns.( |5~6D , (|5.25|) 



f(t, x) = /(to, x) + \f dt' f dydy'a(t', y) dD W>*> x ) y >) 

2 Jt Jf ox 

+ f ' dt'a(t',x)[r ' dy'D(x,y')a i {t',y') + A(t',x)} (5.28) 

Jt Jt' 

6 Futures and Option Pricing 

We derive the futures and options pricing using quantum field theory. For 
the Two-Factor model all the expressions can be obtained exactly; the results 
for the fj, = limit are the same as the HJM-model; we will explicitly give 
the results only for the case of ll — » oo because the expressions for general \i 
don't add much to ones understanding. 

Equation ( |3.7|) for the futures price T only changes for Q; from eqn. (|3.8|) 
we have (note different domains TZ and T below) 



- f DAe -J^dxa(t, X )A(t,x) e J T C[A] ^ 

Z J 

1 ft* f T 

exp{- / dt I dxdx'a(t,x)D(x,x';t,T)cr(t,x)} (6.2) 

2 Jtn Jt, 



and, using the no arbitrage condition ( p. 25 ) we obtain the generalization 
of (|3.10|) given by 



£V(io, t*, T) = — \ dt dxai(t,x) / dx'D(x,x';t,T)ai(t,x') (6.3) 
Jt Jt Jt, 
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For the One-Factor Model with only o\ ^ we have 



lim n^(t ,U,T) = d^-T)(2T-t -Q+O(^) (6.4) 

4/i fl 2 

For the price of a European call option C, a calculation similar to the one 
carried out in Section 4 gives the same formula for *&(G) given in eqn.( |4.14] ) 
with q 2 given in eqn. ([4.12|) replaced by 

rU rT 

q 2 = dt dxdx'(Ti(t,x)D(x,x';t,T)<Ti(t,x') (6.5) 

We have 

rU rT 

jimg 2 = ^ dt(T-t)J dxa 2 (x,t) (6.6) 
For the Two-Factor model we have 



er 2 

lim q 2 = -L(T - U)(U - t ){2T - t - U 



+ 



^|[{l + 2A(T-t )}(l- e - 2 ^-^) 
2A(t,-t )(l-e- 2A(T -*- ) )l (6.7) 



Note for both the futures and option prices, the presence of \i is like adding 
another factor to the model. However, the dependence of the derivatives on 
fi is quite different from that on A; for instance the no arbitrage condition 
changes significantly as \x goes from small to large whereas no such effect 
happens in the case of A; the prices of the derivatives also show non-trivial 
dependence on //. 

If we are interested in pricing any path dependent option or other deriva- 
tives, it is not sufficient to know only the propagator D(x,x';t,T); the full 
structure of the action S[A] is then required. 

For example the payoff function of an Asian option at time to on a zero- 
coupon bond P(t, T) with exercise time £* is given by 

g[P(*, T)} = (— L- ! U dtP(t, T) - K) + (6.8) 

C* — to Jt o 
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Another example is the price of a European call option on a coupon bond 
B(t,T) given in (|2.1| ); the payoff function is then 



g[B] = (B(U,T)-K) + (6.9) 

The payoff function g[A] in both the cases above is path dependent. Ex- 
pressing all the zero-coupon bonds in terms of the quantum field A(t, x), the 
prices of such path dependent options at time to are given by 

C(t ,U,T,X) = ^J DAe-^ dtr(t) g[A]e s ^ (6.10) 

The computation above can only be performed numerically |2(J; for this 
the functional integral over A(t, x) has to be discretized, and which is briefly 
discussed in Appendix A. 

7 Conclusions 

We have re-formulated the theory of Treasury bonds in terms of path integra- 
tion. The HJM-model has a simple path integral with an ultra-local action. 
The statements about martingale conditions and the evaluation of futures 
and options were shown to be calculable in a straightforward manner using 
path integration. The motivation for re-deriving the well-known results of 
the HJM-model was firstly to understand the path integral formulation of 
the quantities of interest in finance, and secondly, to then generalize these 
quantities to the case of quantum field theory. 

The quantum field theory of Treasury bonds is more general than the 
HJM-model; in particular, the correlation of fluctuations of the forward rates 
can be easily modeled to be finite in the field theory whereas in the HJM- 
model all the fluctuations are exactly correlated . From the point of view 
of finance, it is unreasonable to assume that the all forward rates fluctuate 
identically as in the HJM-model; the multi-factors in HJM-model try and 
capture the finite correlation in the time-to-maturity that should exist for 
the forward rates. 

We considered a Gaussian model for the field theory generalization of the 
HJM-model as this is the simplest extension, and also because the model 
could be solved exactly. In particular, the formulae for the futures, cap 
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and option price of Treasury bonds were derived and involved nontrivial 
correlations in the volatility of the model. 

We can generalize the model to account for stochastic volatility of the 
forward rates. This entails introducing another quantum field for modeling 
the fluctuations of volatility, and is similar to the quantum mechanical treat- 
ment of volatility for a single security ||. Stochastic volatility makes the 
system highly nonlinear and is treated in some detail in f23|. 

The best way of modeling Treasury bonds in practice is a computational 



and empirical question [£1], P2| ; only if the field theory model can be easily 
calibrated and yields more efficient algorithms, will it it be taken seriously 
by the practitioners of finance. For the more theoretical side of finance, 
the methodology of field theory certainly adds to the ways of studying and 
understanding the stochastic processes which drive the capital markets. 
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Appendix 

A Lattice Formulation 

We do a more careful and rigorous treatment of the field theory for the 
Treasury bonds. We first discretize the variables into a lattice of discrete 
points. Let (t, x) — > (me,na), where e is an infinitesimal time step and a is 
an infinitesimal in the x direction. Consider the trapezoidal domain T given 
in Fig.(Q) to be bounded by integers m = m (= m = m*(= —), me = na 
and n = N(= -). The integers then take values in the lattice version of the 
trapezoidal domain, say T> given by 

T> = {m = mo, mo + 1, ...m*; n — m, m + 1, ....N} (A.l) 
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The forward rates and quantum field yield on discretization 



f(t,x) -> f(me,na) = f mn ( A - 2 ) 
A(t,x) — > A(me,na) = A mn (A. 3) 



and similarly for a and cr. 
From eqn. (|5.6| ) we have 



fm+ln fmn £&mn ^-®mnA mn (A. 4) 

Using finite differences to discretize derivatives, the generalization of ac- 
tion So in eqn.( |2.10| ) is given, for s = V Na — me , by 

e m * a N la N ~ x 

S[A] = "2 J£ +1 ( (N + l)a - me f/ L + N^^J* J^** 1 " ^ 

(A.5) 



^ = n J o Tl n t/^V r dA*. (a.6) 



ra=mo+l 



Note the functional integral over the field A(t, x) has been reduced to a 
finite- dimensional multiple integral over the A mn variables, which in the case 
above consists of (m* — m ){N — (m + m* — l)/2} independent variables; 
hence all the techniques useful for evaluating finite dimensional integrals can 
be used for performing the integration over A mn . 

To achieve the correct normalization, one in fact need not keep track of 
all the tedious pre-factors in ( |A.6| ). Instead one simply redefines the action 
by 



e s[A] ^ e s[A ]/z (AJ) 
Z = J alAe 3 ^ (A.8) 

All the pre-factors in ( |A.6| ) cancel out; and more importantly, the expres- 
sion e s /Z is correctly normalized to be interpreted as a probability distribu- 
tion, and hence can be used for Monte Carlo studies of this theory. The action 
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given in (|A.5 ) is the starting point for any simulations that are required of 



the model including the pricing of path dependent derivatives; there are well 
known numerical algorithms developed in physics for numerically studying 
quantum fields |20| . 



We explicitly solve for the case of /i — > to see how the HJM-model 
emerges. For fi — > 0, the second term in the action gives a product of 5- 
functions and we have 



e 



m* N-l 

s[A] = e s -q Yl S(A mn+1 - A mn ) (A.9) 

77i=mo+l n=m 

s " = -\ t (N+ ; a _ mf j t<« (a.io) 



Consider evaluating a typical expression like Z in Q5.14 ). For each m, 



there are N-m+1 integration variables A mn ; from eqn.( |A.9| ) we see that there 
are N — m 5-functions, leaving only one variable, say A mm unrestricted. For 
simplicity, we take e = a; hence we have 



m * l~T~ r 

z = II Jir dAmm€ ( A - n ) 

771=7710 + 1 ^7T J 

So = — E A lm (A- 12) 

771=7710 + 1 



Defining W(m) = A mm , we see from eqns. ( p.9| ) that we have recovered 
the HJM-model. We can equivalently consider 



i 



N - m + 
and we have 

\im W(m) -> A mm (A. 14) 

fi— »o 

Taking the continuum limit, we see that the field theory, in the /i — > 
limit reduces to 
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S -> -\ f** dtW 2 {t) (A. 15) 



2 



to 



W(f) = tJ— f T dxA(t,x) (A. 16) 



T-t 

For the general case of fi ^ 0, from eqn.(|A.5|), taking the continuum limit 



of e — > 0, a — > we finally obtain 

/ DA = J] /<M = lim JJ f dA mn (A.18) 
DAe s[A] (A. 19) 



B Generating Functional Z[J] 

Since the generating functional Z[J] has been of central importance in study- 
ing the quantum field theory, for completeness we briefly discuss its deriva- 
tion; all these results are well-known in physics and this derivation is 
intended for readers from other disciplines. 
Recall 



Z[J] = - I DAe s[A ' J] (B.l) 
Z J 

S[A,J] = /'* dt [ T ' dxJ(t,x)A(t,x) + S[A] (B.2) 

Since S[A, J] is quadratic functional of the field A(t,x), to perform the 
functional integration over the field, all we need to do is to find the specific 
configuration of A(t,x), say a(t,x) which maximizes S[A, J]; due to our 
choice of normalization Z[J] depends only on a(t, x). 

Since there is no coupling in the time direction t, we study the solution 
a(t, x) separately for each t, and on the finite line interval t < x < T. We first 
study the case for which the boundary values of the field A(t, x) are fixed, 
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that is consider A(t, t) = p and A(t,T) = p' to be held fixed; we will later 
integrate over p,p' as is required for the evaluation of Z[J\. We henceforth 
suppress the time variable t for notational convenience. 

The 'classical' (deterministic) field configuration a(t, x) = a(x) is defined 

by 



lA{t~t) ~ ° (B - 3) 
a(x = t) = p; a{x = T)=p' (B.4) 

Doing a change of variables A(t,x) = B(t,x) + a(t,x) and a functional 
Taylors expansion we have, from eqn.( |B.3| ) 



S[A + a,J] = S d [a, J] + S[B] (B.5) 

where due to boundary conditions given in eqn.( [B.4| ) S[B] is independent 
of p,p',J. The functional integral over the B(t,x) variables gives only an 
overall constant which we can ignore and hence we have 

1 r+oo 

Z[J] = - / dpdp'e ScM (B.6) 



z 

We now determine a(x); from (|B.3|) we have 



' () '" {r) -a(x) + (T-t)J(x) = (B.7) 



/i 2 dx 2 

a(t) = p,a{T) = p'-t < x <T (B.8) 

Since eqn.( |B.7|) is a linear, the solution for a(x) is given by a sum of 
the solutions of the homogeneous and inhomogeneous equations; it can be 
verified that, for (3 = T — t, we have 



a(x) = . f [a H (x) + a IH (x)] (B.9) 

(B.10) 
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with the homogeneous solution given by 



Q>h(x) = psinhfj,(T — x) + p'sinh/i(x — t) (B-H) 
and the inhomogeneous solution given by 

[T 

o>ih(x) — AW dx'[9(x — x')sinh/j,(T — x)sinh/j,(x' — t) 

+ 9(x' - x)sinhn(T - x')sinhfi(x -t)]J(x') (B.12) 

The 'classical' action is given by 

S cl [a,J] = S 1 [p,p';J] + S 2 [J] (B.13) 

with 

s = - ttfJnhW {-W)(p 2 +/) - 2pp>} + -^y^ 

(B.14) 

/■T /-T 

P = dxsinhfx(x — t)J(x),Q — / dxsinhfi(T — x)J(x) (B.15) 
Jt Jt 

and 



^Pl — T7 — m / dxdx'6(x — x')sinhjj,(T — x)sinhjj,(x' — t)J(x)J(x') 

sinh(fi i) Jt 

(B.16) 

Performing the Gaussian integrations over p,p' and restoring the time 
variable t yields 



Z[J\ = ^e s ' 2 ^ J dpdp'e Sl[p ' p '-> J] (B.17) 

= exp^ /** dt [ T dxdx'J(t,x)D(x,x';t,T)J(t,x f ) (B.18) 
2 Jto Jt 
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where, from eqns.( [B7T3| ), ( |B.14| ) and ( |B.16| ) we have 



D(x,x';t,T) 



sinh d (/i/3) 



sinh ji(T — x) sinh — t){l + sinh 2 (/i/?)6 l (x — x')} 



+ sinh /i(T — x') sinh fx(x — t){l + sinh 2 (/i/5)6'(x / — x)} 
+ cosh(yU/3) {sinh p,(x — t) sinh — t) 
+ sinh /i(T — x) sinh /x(T — x')} 



(B.19) 
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